The Gravitomagnetic Influence on Gyroscopes and on the Lunar Orbit 
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Gravitomagnetism — a motional coupling of matter analogous to the Lorentz force in 
electromagnetism — has observable consequences for any scenario involving differing mass currents. 
Examples include gyroscopes located near a rotating massive body, and the interaction of two or- 
biting bodies. In the former case, the resulting precession of the gyroscope is often called "frame 
dragging," and is the principal measurement sought by the Gravity Probe-B experiment. The lat- 
ter case is realized in the earth-moon system, and the effect has in fact been confirmed via lunar 
laser ranging (LLR) to approximately 0.1% accuracy — better than the anticipated accuracy of the 
Gravity-Probe-B result. This paper shows the connnection between these seemingly disparate phe- 
nomena by employing the same gravitomagnetic term in the equation of motion to obtain both 
gyroscopic precession and modification of the lunar orbit. Since lunar ranging currently provides 
a part in a thousand fit to the gravitomagnetic contributions to the lunar orbit, this feature of 
post-Newtonian gravity is not adjustable to fit any anomalous result beyond the 0.1% level from 
Gravity Probe-B without disturbing the existing fit of theory to the 36 years of LLR data. 
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Part of the post-Newtonian gravitational interaction 
between two mass elements, when both are in motion, 
has been called "gravitomagnetism," in analogy with the 
magnetic force between moving charges. The gravito- 
magnetic interaction is part of the more general 1/c 2 - 
order motional corrections to Newtonian gravity that re- 
sult from field theories such as Einstein's general rela- 
tivity and scalar-tensor generalizations P, 13| ■ A total 
package of velocity-dependent corrections is required so 
that the gravitational equation of motion remains consis- 
tent when expressed in different asymptotic inertial ref- 
erence frames. If Lorentz invariance of local gravitational 
physics is imposed by empirical constraint, the package of 
motional corrections is additionally limited in structure. 

The gravitomagnetic interaction of general relativity 
was first studied by Lense and Thirring in 1918, and it 
was shown to produce both accelerations of and torques 
on two neighboring rotating bodies. Others, viewing this 
phenomenon geometrically, have coined the interpretive 
name "inertial frame dragging" from rotating matter. It 
has also been shown that the gravitomagnetic interaction 
plays a part in both shaping the lunar orbit at a level 
art in a thousand) readily observable by laser ranging 
, and in contributing to the periastron precession of 
binary and especially double pulsars [J] . 1 



1 Although the present analysis uses the equation of motion from 
the parameterized post-newtonian (PPN) formulation of long 
range, metric gravity, the conclusions we draw are not con- 
strained by this choice. The phenomenology we explore — the 
Vi X (vj X gij) acceleration of Eq. JTJ — generically results from 
modified metric field expansions or even non-metric models of 
gravity at the 1/c 2 post-Newtonian level. If such an interaction 
is used to explain precessional effects in a gyroscope experiment, 



For applications to the analysis of gravitational phe- 
nomena, a general metric tensor field expansion for the 
gravitational potentials in a broad class of theories was 
developed by Will and Nordtvedt 0, 0] . This parameter- 
ized post-Newtonian (PPN) framework yields a gravito- 
magnetic contribution to the equation of motion, which 
in the Lorentz-invariant case is 

a t = (2 + 27)^^3- Vi x ( Vj x ry), (1) 
3 C r « 

where Vi and Vj are the velocities of bodies i and j in 
the chosen asymptotic inertial coordinate system Q . The 
vector rij, when combined with the fraction [tj/rfj con- 
stitutes the Newtonian gravitational acceleration of mass 
i toward mass j. In geometric language, the PPN fac- 
tor 7 quantifies the amount of space curvature produced 
per unit mass. In general relativity, 7=1. Metric the- 
ories allowing preferred inertial frame effects (absence of 
local Lorentz- invariance) add the parameter ai/4 to the 
(2 + 27) pre-factor in Eq. |T]), but lunar laser ranging as 
well as other solar system observations constrain a\ to 
be less than 10" 4 @]. 



then it will also be present to perturb the lunar orbit and bi- 
nary pulsar orbits, regardless of its parameterized strength in any 
particular model. Any attempt to suppress the gravitomagnetic 
influence on the lunar orbit relative to that on a low-orbiting 
gyroscope by a Yukawa-like modification to — whether met- 
ric or non-metric — would clash with strong constraints on the 
inverse-square nature of gy determined via satellite and lunar 
laser ranging. Eq. JTJ is, of course, dependent on the asymptotic 
inertial frame in which analysis is performed, just as are magnetic 
forces within an electromagnetic system. Consistent formalisms 
can be used in any choice of frame to calculate observables; this 
permits using the convenient solar system barycentric frame for 
our analysis. 
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We can ask what effect the gravitomagnetic term of 
Eq. (1} has on a gyroscope outside of a rotating spherical 
mass. We define the gravitomagnetic field by 



G ij = (2 + 2 7 )J2 



c 2 r 3 



(t>j X ry 



(2) 



so that a,i = Vi x dj in analogy to the electromag- 
netic Lorentz force. Considering a small gyroscope, the 
Gij vector field is calculated at the gyroscope center, 
and will be nearly constant across its body. To obtain 
the cumulative effect of mass elements moving within a 
body rotating at angular velocity Q, the gravitomagnetic 
field is integrated over all mass elements, j, each with 
d/ij = Gp(rj)d 3 rj, where G is Newton's gravitational 
constant, and p(rj) is the mass density at radius Tj from 
the body center. Adopting a spherical coordinate system 
aligned with the rotation axis of the body, we describe 
the Cartesian vector rj = r sin 9 cos <j>\ + r sin 9 sin 0j + 
r cos 9k, and the vector to the gyroscope (placed in 
the = plane) is = asin^i + acosipk, so that 
rfj = a 2 + r 2 — 2ar (sin tp sin 9 cos + cos ip cos 9) . The ve- 
locity of mass element j is Vj = fir sin 9(— sin 0i+cos 0j) , 
so that 



dGa = (2 + 2 7 ) 



fir sin 9d[ij 
c 2 r 3 - 



x {(r cos9 — a cos t/0 ( cos + sin0j) 

+ (asini/>cos0 — rsin0)k}. (3) 

For the special case of a gyroscope situated over the 
pole for simplified integration (ip = 0), and recognizing 
that the k component of the Gij vector will be the only 
one to yield a non-zero angular integral, we find that 



(2 + 2 7 )Gfik 



l Jo (a 2 + r 2 — 2aru) ? 
s, ( 2-;-2 7) ^k r R drp{r y, ( 4 ) 



3c 2 a 3 



Here we used the identity u — cos^, and note that the 
integral over u eliminates the r-dependence in the de- 
nominator. Recognizing that the moment of inertia of a 
spherical body is 



we see that 



GiM = 0) 



(2 + 2 7 )G/ s n 



(5) 



(6) 



where the s-subscript represents the massive rotating 
body. 



Orienting the gyroscope so that its spin axis is along 
the i direction: oj — io\, the velocity of an element within 
the gyroscope is =wxr, = — cor cos 9] +ur sin 9 sin 0k, 
where is the vector position of a mass element within 
the gyroscope with respect to its center. Now all pieces 
are at hand to evaluate the acceleration of each mass ele- 
ment within the gyroscope due to the gravitomagnetic 
term. The force on each element is dFi — dwiidi = 
drriiVi x Gij, and the torque on the gyroscope from this 
element is then g?t = x dFi. Combining these steps, 



dr = 



(2 + 2"f)GI s flwr 2 cos 9dm l 



(cos 9j — sin 9 sin 0k) . 



(7) 

Integrating this torque over the volume of the gyro- 
scope, the k component integrates to zero in the inte- 
gral, yielding 



(2 + 2j)GI s nuj . 
T = 37s J 



(l+ 7 )GW g w, 



1 r-R 

u 2 du I p(r)r 4 dr 
-l Jo 



(8) 



where the integral is seen to be one-half the rotational 
inertia of the gyroscope, which is denoted with subscript, 
g. This torque will change the angular momentum vector 
of the gyroscope, L g = I g uj, such that the angle of the 
axis, $, precesses at a rate of $ = |t| / \L S \, so that 



(1 + 7 )GI s fi 



(9) 



The direction of precession indicated by Eq. (JSJ) is 
one that takes the angular momentum — originally in the 
i direction only — toward the j direction, meaning that 
the precession has the same sense as the rotation of 
the massive body. Had we developed an expression for 
G^ (ip = f ) at the equator of the rotating body, we would 
have found half the magnitude of the polar case and in 
the opposite direction. In general, the field 



G, 



(i + j)Gi s n 



3(k • r)f - k] 



(10) 



where f is the unit radial vector. In a circular polar orbit 
with uniform angular rate in tj), the gravitomagnetic field 
averages to 



{Gij) 



(i + 7 )Gj s n 

2c 2 a 3 



k, 



(11) 



so that the net field shares the same direction as that 
over the pole, and therefore the net precession will be in 
the same sense as the rotating mass, but at one-fourth 
the polar rate. Summarizing, 



polar orbit 



(l + 7 )G/ S rj 

4c 2 a 3 



(12) 
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Putting this in another form, where we reduce the ro- 
tational inertia to I s = fM s R 2 , where / = 0.33 for the 
earth [8J], we have the more convenient form: 



polar orbit 



(1 + 7)/ GM. 
4 Rc 2 



a 



(13) 



For Gravity Probe-B (GP-B), in a 640 km altitude polar 
orbit, Eq. (fT3"|) yields 0.042 arcseconds per year, matching 
the mission expectation 0, [l(| • GP-B anticipates mea- 
suring this precession to < 1% accuracy, and perhaps as 
well as 0.1-0.3%. Note that the gyroscope spin was not 
treated as an intrinsic property in deriving the gyroscope 
precession. Rather, the effect results from the integrated 
mass currents of mass elements in rotational motion. 

In obtaining the effect of the gravitomagnetic term 
(Eq. ((TJ) on the lunar orbit relative to Earth, we treat the 
gravitomagnetic acceleration as a perturbation about an 
otherwise circular orbit. We start with an orbit obeying 
the two-body equation of motion: 



GM 



o(r) 



(14) 



where a(r) is the central acceleration, and I = r 2 u is 
the (specific) angular momentum. We idealize the un- 
perturbed orbit to have zero eccentricity and zero incli- 
nation, so that the end result is accurate for the moon 
at the 10% level or better. The deviation, Sr, resulting 
from a periodic acceleration perturbation, Sa, then obeys 



Sr 



Sa r + 2cjq / Sa T dt' 



(15) 



where ujq is the natural frequency for orbital perturba- 
tions, with ujq w 3ui 2 — jr ~ uj2 - The acceleration, Sa, has 
been decomposed into radial and tangential components, 
and t' is a time variable. 

Expressing the triple cross-product in Eq. (JTJ as the 
equivalent dot-product relationship, we find that the 
gravitomagnetic acceleration of the moon is 

a m = Sa = ^ e ( 2 + 2 T) [f( Vjn . Ve ) - Ve (v m ■ f)} , (16) 
c z a z 

where r is the unit vector from the earth to the moon, 
and a is the earth-moon distance. Eq. p6[) is re-written 
as 



V ■«) -V(V -f + u-f)] , 
(17) 

with the earth's velocity around the sun being V and the 
moon's velocity being V + u, where u is approximately 
thirty times smaller than V in magnitude. 

Note that u represents to sufficient accuracy the 
geocentrically- viewed orbital velocity of the moon around 
the earth. Thus under the assumption of a circular orbit 
(about which we examine the perturbation), u • r = 0. 



Likewise, if we define r to be the tangential orbit vector 
at the moon that is perpendicular to r, u ■ t = u. Under 
the assumption that the earth is in a circular orbit about 
the sun, the relationship between V (perpendicular to 
earth-sun line) and r and f picks out the synodic phase 
angle of the moon, D. Specifically, V ■ r = —V sin D and 
V ■ t = — VcosD. Similarly, V • u = —VucosD. Leav- 
ing off the pre-factor from Eq. (117|) for now, and dealing 
only with the vector math, the radial component is 



Sa r cx V 2 — VucosD — (- 



-V 2 cos2L> 
2 



and the tangential component is 



VsinD) 2 
- VucosD, 

1 



(18) 



Sa T oc -(-VcosD) (-VsinD) = --V z sin2D. (19) 

The periodic accelerations consist of two categories: V 2 
terms that have 2D angular dependence, and Vu terms 
that have D angular dependence. We can treat each 
separately in solving Eq. (Tl3|) . There is also a constant 
term in the expression for Sa r . We can ignore the con- 
stant term since it only acts to rescale the orbit in a 
non-periodic way. We will deal first with the 2D terms, 
then look at the D terms. 

First, we integrate the Sa T term. Noting that the rate 
at which D advances is D = uj — Cl, the difference be- 
tween the lunar orbital frequency and the earth's orbital 
frequency, we can construct an arbitrary 2D argument as 
[2(uj — Q)t' + </>], where <j) is an arbitrary phase depend- 
ing on the choice of t' = 0. The integral (without the 
numerical pre-factor) is then 



2uj / sin[2(w - Q)t' + <p]dt' 



-n 



cos 2D + const. 

(20) 

We have consolidated any initial phase in the integration 
constant, effectively defining t so that D = (lo — VL)t. As 
above, we can ignore the constant term in our periodic 
analysis. The differential equation becomes 



Sr + uj Sr 



(l +1 )GM v2 



cos 2D 



- fl 



■ cos 2D 



(1 + 7)GAf V 2 2uj-n 



n 



cos 2D. 



(21) 



The solution, in meters, is then 
V 2 2 - r) a 



Sr w -(l+ 7 )^_ 1- — _cos2£> » -6.5 cos 2D m, 

c z 1 — 77 3 — 877 

(22) 

where we have made use of Kepler's relation (uj 2 a 3 — 
GM) and define 77 = Q/ui, ignoring terms to second order 
in n. 

The term proportional to Vu has no tangential part, 
so we immediately write the differential equation as 



Sr + uj Sr 



(2 + 2 7 )GM 



Vu cosD, 



(23) 
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for which the solution is 

(2 + 2 7 )GM Vu 



5r = — - 



c 2 a 2 to 2 - (w - n) 2 



cos D 



-(1 + 7) 



Vu u 



a cos D w -3.4 cos D m. (24) 



But a feedback process produced by the interaction of 
synodic perturbations and the cos 2D solar tidal distor- 
tion of the lunar orbit results in an amplification of cos D 
terms by the factor 11 1 



1 - 2r\ 



1.79, 



— 1 - 7 V 
so that the corrected range oscillation is 



(25) 



Vu to 1 — 2i] 

Or ss — (1 + 7)^5--=- a cos D w —6.1 cosD m. (26) 

&■ s2 1 — Y77 

In summary, the gravitomagnetic perturbations of the 
lunar orbit are: 



8r2D = —6.5 cos 2D meters 
5ro = —6.1 cos D meters. 



(27) 



Lunar laser ranging (LLR) has been used for decades 
to provide a number of the most precise tests of general 
relativity, including tests of the weak and strong equiv- 
alence principles, time-rate-of-change of Newton's grav- 
itational constant, G, geodetic precession, among oth- 
ers Equivalence principle violations would produce 
a cosD signal 13], though no cos 2D signal. Current 
fits to the archive of LLR data limit any net devia- 
tion of the cos I? term in the lunar orbit to less than 

4 mm from the orbit prescribed by general relativity 
[l2| . Likewise, the cos 2D term is constrained at roughly 
the 8 mm level. Thus barring a chance simultaneous vio- 
lation of the equivalence principle and gravitomagnetism, 
the 4 mm constraint translates to a check on the ^6 m 
gravitomagnetism amplitude to better than 0.1% accu- 
racy. Allowing for such a conspiracy, we must use the 
8 mm cos 2D constraint (which is not influenced by equiv- 
alence principle violation) to establish a » 0.15% verifi- 
cation of the gravitomagnetic phenomenon. At this time, 
LLR provides the most precise test of this phenomenon — 
far better than the LAGEOS tests of the Lense-Thirring 
effect [13] and tests from binary pulsars Q . This result is 
also likely better than the ultimate result from the GP-B 
experiment [9|]. 

A new effort in LLR is poised to deliver order-of- 
magnitude improvements in range precision [po] ]. which 
will translate into tighter constraints on the cosD and 
cos 2D amplitudes in the lunar orbit. Because these 
are periodic effects, their accurate determination requires 
only about a year of new data collection. Thus a signifi- 
cantly improved test of this phenomenon is not far away. 



Whether the mass elements in a body are moving 
commonly — as for the Earth and the Moon in orbital 
motion in the solar system — or as mass currents in the 
rotational manner of the spinning Earth and gyroscope 
in GP-B, the total interaction between bodies must be 
dominated by a linear-order integration over the bod- 
ies' mass elements in both situations. Breaking weak- 
field superposition would be a radical and ultimately 
non-viable modification to gravity theory. If this linear- 
order gravitomagnetic interaction from PPN metric grav- 
ity, Eq. ([T]) , is altered in order to fit any anomalous GP-B 
observation, then either the cosD or cos 2D amplitudes, 
or both, in Earth-Moon range as measured by LLR to 
half-centimeter accuracy will show anomalies under this 
new modeling — establishing a profound empirical clash. 
An added likely consequence of modifying gravitomag- 
netism would be destroying the total fit to the binary 
pulsar 1913+16 data which includes a better than one 
percent match to General Relativity's predicted gravi- 
tational radiation-reaction accelerations in that system. 
Existing and robust observations already encumber the 
gravitomagnetic interaction. 
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